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Abstract 

We study numerically the mechanical stability and elasticity properties 
of duplex DNA molecules within the frame of a network model incorpo- 
rating microscopic degrees of freedom related with the arrangement of 
the base pairs. We pay special attention to the opening-closing dynamics 
of double-stranded DNA molecules which are forced into non-equilibrium 
conformations. Mechanical stress imposed at one terminal end of the DNA 
molecule brings it into a partially opened configuration. We examine the 
subsequent relaxation dynamics connected with energy exchange processes 
between the various degrees of freedom and structural rearrangements 
leading to complete recombination to the double-stranded conformation. 
The similarities and differences between the relaxation dynamics for a pla- 
nar ladder-like DNA molecule and a twisted one are discussed in detail. 
In this way we show that the attainment of a quasi-equilibrium regime 
proceeds faster in the case of the twisted DNA form than for its thus less 
flexible ladder counterpart. Furthermore we find that the velocity of the 
complete recombination of the DNA molecule is lower than the velocity 
imposed by the forcing unit which is in compliance with the experimental 
observations for the opening-closing cycle of DNA molecules. 
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1 Introduction 



Nowadays powerful micromanipulation techniques allow for the selective exci- 
tation of DNA in single molecule experiments under controlled conditions (see, 
e -g-> D>E3>E9)- Particularly mechanical properties of DNA molecules under 
tension and the inducement of structural transitions with external stress have 
received a lot of attention In this context several force measurements 

were performed on single DNA molecules to examine their elastic response to 
applied stretching and twisting forces, respectively 0] , , [H] . It has been 
demonstrated that with the application of stronger forces to the molecule in 
stretching experiments, even highly cooperative conformational transitions were 
discovered, for which, e.g., B-DNA is converted into a new overstretched con- 
formation called S-DNA 6 . Furthermore, with sufficiently high tension the 
mechanical unzipping of single DNA molecules is achievable jo]. The opening 
process of DNA has been the subject of a number of studies due to its immense 
relevance for the transcription process. Simple models taking into account the 
elasticity properties of DNA, such as the wormlike and rodlike chain model, 
respectively have been invoked to describe these unzipping experiments |1U| . 
The statistical mechanics of the torque induced denaturation of DNA has been 
studied in Recently, a semi- microscopic model of the binding between the 
two strands has been used to explain the effects of applied tension to DNA 
unzipping JT2] . In order to determine the geometry and deformability of DNA 
at a local level, i.e. in terms of base pairs, energy functions based on discrete 
network models have been derived in |14 j -jl7 j . In an early network model for the 
base pair opening in DNA involving microscopic degrees of freedom, Peyrard 
and Bishop (PB) proposed a planar ladder-like oscillator model of DNA assign- 
ing each base pair a vertical inter-strand vibrational degree of freedom which 
simulates the stretchings and compressings of the corresponding H-bridges |18|. 
The binding forces of the hydrogen bridges are described by a Morse poten- 
tial. The bases itself are treated as point masses. Horizontally, the bases on 
the same strand are coupled via the stacking interaction described by harmonic 
potentials. The PB model has been extensively studied and localized oscillating 
solutions (breathers) have been found reflecting successfully some typical fea- 
tures of the DNA opening dynamics such as the magnitude of the amplitudes 
and the time scale of the breathing of the 'bubble' occurring prior to thermal de- 
naturation O n the other hand, as the coiled form of double-stranded DNA 
is concerned, it has been found that the bubble formation is strongly correlated 
with twist deformations and local openings are always connected with a local 
untwist of the double helix , [201 • I n order to account for the hclicoidal struc- 
ture of DNA the PB model has been significantly extended by Barbi, Cocco and 
Peyrard (BCP) 19 . In their model of the twisted form of DNA (the Watson- 
Crick double helix), two degrees of freedom per base pair are introduced. There 
is a radial variable measuring the distance between two H-bridged bases along a 
line that connects them in the base plane being perpendicular to the helix axis. 
Further, the twist angle between this connecting line and a reference direction 
determines the orientation of the H-bridge. In this model a further degree of 
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freedom has been included, namely the axial distance variation Interac- 
tions between the bases are described then by appropriate potential terms in the 
energy of the system. The process of DNA unzipping by applied force has been 
studied in great detail ^3,^21iC31- The merit of these works is that many of 
the theoretical conclusions, gained on the basis of rather simple network models, 
are in agreement with the experimental results. 

In the current study of the nonlinear dynamics related with the opening pro- 
cess and elasticity features of DNA molecules our aim is to examine the energy 
exchange processes and the relaxation dynamics in DNA molecules after they 
have been brought into a non-equilibrium conformation by applied tension. Such 
an investigation is associated with recent mechanical experiments performed 
with sing le DNA molecules 0-0 forced away from their equilibrium confor- 
mations. Following force applications energy redistribution within the DNA 
molecule takes place such that a (new) equilibrium conformation is attained 
[21] [23 We focus our interest on the elasticity properties and the relax- 
ation dynamics of two different double-stranded DNA configurations, namely 
the planar untwisted (ladder-like) system and the twisted one arising when the 
two strands of the molecule are coiled around its molecular axis. In particu- 
lar we address the questions in which way these two different DNA molecules 
react to imposed stress and how fast the equilibration process takes place for 
them. Furthermore, we compare the original equilibrium structures with those 
newly attained after release from the exerted mechanical stress. Regarding the 
biological functioning of DNA their elasticity properties are of importance for 
instance for the creation of a denatured bubble appearing prior to strand separa- 
tion making the DNA accessible to the transcription and replication processes. 
Compared to previous theoretical studies of the unzipping dynamics of DNA 
induced by strong forces yielding complete strand separation JJ,^21 we re ~ 
mark that for the current study we rather concentrate on mechanically stressed 
double-stranded DNA which gets forced away from the equilibrium configura- 
tion, however, only with beginning strand separation. From a biophysical point 
of view this also resembles the situation when the DNA molecule is inflicted 
to internal stress changes occurring in the living cell which alone is not strong 
enough to cause advanced bond breaking. Furthermore, the interaction be- 
tween DNA and proteins such as enzymes can exert forces that initiate the local 
opening of the DNA molecule. 

The paper is organized as follows: In the second section we describe our 
network model for the structure of double-stranded DNA. The third section 
deals with the relaxation dynamics within DNA molecules forced into locally 
distorted configurations. We study the similarities and differences between the 
relaxation dynamics for a planar ladder-like DNA molecule and a twisted one. 
Finally, we summarize our results. 
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2 Network model for the double-stranded DNA 

Our DNA model takes the basic geometrical features of the double-stranded 
form of DNA, as a polymeric molecule composed of two strands of nucleotides, 
into account which are essential to model structural deformations generated by 
mechanical stress imposition and vibrational relaxation processes. Likewise the 
approach in ^2>E3> we treat the double-stranded DNA as a network of cou- 
pled oscillators incorporating essential microscopic degrees of freedom of DNA 
and the inherent interactions between them. Such a network model approach is 
appealing because it can successfully describe the mechanical behavior on the 
microscopic scale relevant for biomolecular processes, but bears enough simplic- 
ity to be tractable with not too much computational efforts. The constituents 
of the oscillator network model represent the nucleotides which are regarded 
as single nondeformable entities. Thus, no inner dynamical degrees of freedom 
of the nucleotides are taken into account which is justified by the time scale 
separation between the small-amplitude and fast vibrational motions of the in- 
dividual atoms and the slower and relatively large-amplitude motions of the 
atom groups constituting the nucleotides [21 ■ Concerning the structural com- 
ponents, each nucleotide is composed of a sugar, a phosphate and a base. The 
sugar-phosphate groups of neighboring nucleotides on the same strand are linked 
via covalent bonds establishing the rigid backbone to the strand. There is a base 
attached to every sugar. Since, for simplicity, we do not distinguish between the 
four different types of bases, the nucleotides are considered as identical objects 
of fixed mass. Two bases on opposite strands are linked via hydrogen bonds 
holding the two strands of DNA together. 

For computational convenience we cast the double helix structure in a Carte- 
sian coordinate system whose z— axis coincides with the central helix axis as 
sketched in Fig. 1. The base pairs are situated in planes perpendicular to the 
central helix axis and the vertical distance between two consecutive planes is 
given by h. For the equilibrium configuration each base has equilibrium coor- 
dinates xf\, and z®\. The index pair (n,i) labels the n-th base on the 
i—th strand with i = 1, 2 and 1 ^ n ^ N , where N is the number of base pairs 
considered. The equilibrium distance between two bases within a base pair, do, 
is determined by 



where d n = x { n [ — x { n ' 2 and d v n — y^ [ — y^ ' 2 are the projections of the line 
connecting the two bases on the x, y— axes of the coordinate system. The devi- 
ations d n from do through displacements, x n ,i, Vn,i an d z Ul i, of the bases from 
their equilibrium positions, x^ ■ , y^ \ and z^\ , are expressed as 



d n — V (d% + x„,i - x n ,i) + (d n + y n ,i - 2M,2) +(271,1-271,2) -do. (2) 



The relative twist between two consecutive base pairs, being rotated around the 
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Figure 1: Left: Geometry of the twisted double-stranded DNA model showing 
the variables. Right: The DNA molecule at the end of the pulling process with 
exaggerated distortions. 



central axis, is given by an angle 8q. For later use we introduce the quantity 

6 n = arctan d } + Vn ' 1 ~ Vn ' 2 + 2 TO7r , (3) 

d x 4- T i — T o 

as the angle between the a;— axis (as the reference direction) and the line con- 
necting two (displaced) bases of a base pair measuring the alignment of the 
associated H-bridge. m is an integer to assure monotonicity of 9 n with respect 
to n. 

The three-dimensional equilibrium distance between two adjacent bases on 
the same strand is given by 



\h\ (4) 



fc = V(*S-A.*) a + (fS-A 

and deviations from Iq are determined by 

ln,i = | ( ^riji x n,i ~ x n— + ( ^n.i 2AM — Vn—l,i) 

+ {h + z n ,i - Zn-\,i) j -lo, (5) 



with L na = X n,i ~ X n-l,i and = Vn\ ~ Vn-l,i 

Our model Hamiltonian is then given by 



H = E km + V h + V c + V s + V l . (6) 



Vh and V c represent the potential energy part for the H-bonds and the covalent 
bonds, respectively. The deformations of the hydrogen bonds corresponding to 
deviations d n (t) from the equilibrium bond length do are described by vibrations 
in a Morse potential which is given by 



N 



u h = D J2 E [ cx p(-" d «) - if 



(7) 



1=1,2 n=l 



D is the depth of the Morse potential and a is the range parameter. For the 
sake of simplicity we do not distinguish between the two different pairings in 
DNA, namely the G-C and the A-T pairs. The former pair involves three 
hydrogen bonds while the latter involves only two. Thus our model applies to 
pure poly(dG-dC) and poly(dA-dT) DNA sequences. 

The covalent bonds between the sugar-phosphate groups of neighboring nu- 
cleotides on a strand are rather strong and rigid (with bond energies of the 
order of 2 — 10 eV) in comparison with the weak and flexible hydrogen bonds 
(with bond energies of the order of 0.04 — 0.3 eV). Therefore, it seems justified 
to treat the potential of the covalent bonds, simulated by elastic rods, in the 
harmonic approximation given by 



JV 



U r = K 



E E z 

1=1,2 n=l 



2 

n,i ' 



(8) 



where K is the elasticity coefficient. 

The potential term V s takes stacking effects into account which impede that 
one base slides over another For the form of V s we adopt the one used in 
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= 77 ^ _ d„_ M ) 2 exp[-/3(d M + dn-i,i) 



(9) 



The decrease of the molecular packing with substantially increased base pair 
opening is accounted for by the exponential factor with attenuation parameter (3. 
As the changes of the axial distance between adjacent base planes is concerned 
we assume that the supposedly small longitudinal helix deformations can be 
modeled by a harmonic elasticity potential term given by 



Vi — — (z n ,i — z n -i,i) z 



(10) 



The kinetic energy is determined by 



E ^ = 2L E E (pS) 2 + W0 2 + (*0 ; 



1=1,2 n=l 



(11) 



where m is the mass of a base and p^f denote the (x, y, z)— component of its 
momentum. 
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In order to assign values to the various parameters we note that the geo- 
metrical parameters of the equilibrium configuration are well known |23| . The 
rotation angle for the twisted configuration is 8 = 36°, the distance between 
base pair planes is h = 3.4 A, and the inter-base distance is do = 20 A. (For 
the twisted configuration do corresponds to the diameter of the helix.) The 
ladder-like system is completely untwisted, viz. 9q = 0. For the average mass 
of one nucleotide we use M = 4.982 x 10~ 25 kg. Like Barbi et al ^5] we set 
a = 4.45 A -1 , D = 0.04 eV, and K = \.§eV k~ 2 . Concerning the parameters 
C and S there is little experimental evidence from which an estimate of them 
could be inferred. In j2H| the value of the parameter S has been set to S = 2K 
which we adopt here. (Alternatively in a model approach of DNA denaturation 
dynamics the parameters S and C were specified to fit the melting tem- 
perature of certain DNA polymers.) The parameter C is treated as adjustable. 
Furthermore we use for the value of the attenuation parameter j3 = 0.5 A -1 

With a suitable time scaling t — » \J 'Da 2 jmt and dividing the original Hamil- 
tonian by D one passes to a dimensionless formulation with quantities: 



Pni 



C 



ax 



a 2 D 



ay. 



UZr, 



Av) 



(x) 

Pn,i ~{y) Pn.i ~( z ) 

i Pn.i / 7T J Pn.i 



, jf-i, s- s 



Pni 



a 2 D 



a 2 D 



d n = a d n , f o = a ro , h = ah , = 







(12) 
(13) 
(14) 
(15) 



In the following, we omit the tildes. 

The equations of motion are derived from the Hamiltonian presented in Eqs. 
iJBJi-ljnj and read as 



A*) 



(y) 



P 



(x) 
Pn.i ' 



= 2 [exp(-dn) - 1] exp(-d n ) 



dd n 



2K 
S 



dlr, 



+ ln+1 



1 ^ 2^ dn " dn ^ 



1 + -^{dn+l - dn) 
dd n 



(16) 



(d n ~ d n -%) exp[-^ (dn + dn-i)] 
(d n +i - d n ) exp[-/3 (d n +i + d n )] 



(17) 
(18) 



= 2 [exp(-d n ) - 1 ] exp(-d n ) 



dd„ 

dy n ,i 
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- IK 

- S 



, dl n ,i , dl n +i t 



dy„ 



dy„,i 



1 - ~~ dn ~^ 



1 + -^{d n +l - &n) 



{d n - d n -i) exp[-/3 (d„ + d„_i)] 
(d n+ i - d„) exp[-/3 (d„+i + d n )]| 



•(«) 



dy n ,i ' 
(*) 

= 2 [exp(-d„) - 1] exp(-d n ) 
- 2K 



(19) 
(20) 



dd n 



81 

I 7 

r hi+l,i 



dz 

dl n +l,i 



- S 



'11,1 rs 

dz n 

1 — ^(^n — d n -l 



1 + ^( d n+l - dn) 

dd n 

C ( 2 Zn+l,i Z n _ l,i ) , 



(92 ■ 

(d„ - d„_i) exp[-/3 (d„ + d n _i)] 
(d„+i - d„) exp[-/3 (d„+i + d n )]| 



with the derivatives 



<9d» 
(9/ 



{-lf +1 {d x n +x n s-x na ) 
d n + d a 

L x -\- X n i X n — l.i 
ln,i + ^0 



(21) 

(22) 
(23) 



and the equivalent expressions for dd n /dy n ,i, dl n ^/8y n s, 8d n /8z n ^ and 8l n s/dz n ^. 

The values of the scaled parameters are given by K = 0.683, ro = 44.50, 
h = 15.13 and l — 31.39. One time unit of the scaled time corresponds to 
0.198ps of the physical time. 



3 Relaxation dynamics of mechanically stressed 
DNA 

3.1 Pulling process in the ladder and twisted configura- 
tion 

We study the energy exchange process taking place after the double-stranded 
DNA has been forced mechanically into a non-equilibrium configuration. In 
more detail, we consider the situation when at one end (referred to as the pulled 
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end) of short DNA molecules consisting of 10 up to 30 base pairs the two strands 
get pulled laterally in opposite directions by an external force (exerted e.g. when 
the two strands of one end of the molecule are attached to two different solid 
supports that are progressively separated [S]) so that they are dragged apart 
leading to gradual opening of the double-stranded chain. The pulled end is 
formed by the base pairs with index n — 1 and the rest of the molecule with 
n > 1 is located above it in the z direction. Therefore a positive value of 
z n ,i — Zn—i,i w ih represent a local ^-stretching of the strand i. In order to 
mimic the pulling process we assume that the pulling force is directed solely 
along the x— direction and the x— coordinates of the terminal bases are shifted 
uniformly with a pulling velocity of v p = 20 fxm/s while the corresponding y and 
z coordinates are kept fixed at = z\j = 0. 

Moreover, due to the constraint imposed by the backbone rigidity there 
result deformations not only at the pulled end but also for the rest of the helix. 
In general, the extent of lateral deformations subsides along the strands from 
the pulled towards the free end. The motion of the pulling unit is assumed to 
terminate when the lateral elongation at the pulled end has reached a certain 
value, taken as d = 2.31 A, so that the corresponding terminal hydrogen bond 
is broken. Moreover, the next-neighboring hydrogen bond is broken too and the 
third one is already considerably stretched. The typical resulting deformation 
pattern of the double stranded DNA is depicted in Figs. 2 and 3 for the ladder 
and twisted configuration, respectively. (In the Figures energy is expressed in 
terms of eV , the unit for length is A and the time is given in ns.) 

Comparing the resulting DNA deformations we note that for both, the ladder 
and the twisted configuration, virtually equal lateral elongations are obtained. 
However, unlike the ladder configuration the twisted one experiences additional 
changes of the twist angles. There arises an overall decrease in the twist angle 
and the further a base pair is apart from the pulled end the more is the corre- 
sponding twist angle reduced (cf. Fig. 3 (b)) leading to an untwist of the helix. 
Nevertheless, with growing distance from the pulled end the difference between 
two consecutive angle changes becomes smaller. Thus the relative twist angle 
between two adjacent base planes diminishes towards the free end. 

While for the twisted system the axial deformations are negligible (see Fig. 3 (c)) 
we note that for the ladder configuration there results noticeable diminution 
of the axial equilibrium distances between neighboring base planes which is 
strongest at the free end and weakens gradually towards the pulled end (see 
Fig. 2(b)). The relations between the partial deformation energies for the twisted 
system are shown in Fig. 4. The major contributions stem from the Morse po- 
tential energy associated with the H-bridge deformations and the stacking in- 
teraction energy both exceeding the energetic content of the covalent bond and 
axial distance change energy, respectively by four orders of magnitude. Typical 
for DNA molecules is that there exist different elasticity regimes, that is the 
degree of the deformation depends not overall linearly on the force (see e.g. [3]) 
which is reflected here in the nonlinear character of the temporal behavior of 
the deformation energies. 



9 



Figure 2: Deformation pattern at the end of the pulling process for the ladder 
system, (a) The lateral deformations d n . (b) The deformations in axial direction 




Figure 3: Deformation pattern at the end of the pulling process for the twisted 
system, (a) The lateral deformations d n .(b) The twist angle changes, <f> n = 
0„ — (n — l)#o -(c) The deformations in axial direction z n . 
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Figure 4: Pulling process for the twisted system: The temporal behavior of the 
partial deformation energies as indicated in the figure. 



3.2 Relaxation processes 

After the DNA molecule has been brought into a non- equilibrium configuration 
with imposed mechanical tension we explore now the vibrational dynamics of the 
helix cage after the pulled end has been released from the forcing lever. Hence, 
there is no external mechanical constraint enforced on the double helix any 
longer. For the then isolated double helix we expect that energy redistribution 
processes between the various vibrational degrees of freedom take place targeting 
for the attainment of a new stationary configuration. 

We discuss the ladder system and the twisted one separately and start with 
the former case for stacking interaction strength S = 1.26. 

We integrated the set of coupled equations l|ltj|) - (|19|) (with free boundary 
conditions) with a fourth-order Runge-Kutta method while the accuracy of the 
computation was checked through the conservation of the total energy. 

To gain insight into the energy redistribution process we monitor the tem- 
poral behavior of the potential energies associated with the deformations of the 
hydrogen and the covalent bonds, respectively, and changes of the axial elasticity 
and stacking interaction energy, respectively. Fig. 5. illustrates how the total de- 
formation energy of 0.35 eV, which has been deposited in the DNA molecule, is 
shared between the various degrees of freedom in the course of time. Generally, 
we observe that in an initial interval of ~ 7000 ns, hereafter referred to as the 
redistribution time, pronounced migration of energy takes place. Regarding the 
Morse potential energy we note that during this initial interval the amplitude 
oscillates wildly in a range reaching from 0.05 eV up to 0.30 eV . Afterwards the 
potential energy performs lower-amplitude oscillations around the mean value 
0.1 eV indicating that a quasi-equilibrium regime is attained. Apparently, the 
amount of energy by which the Morse potential energy gets reduced (almost 
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50 % of its initial value) has been redistributed into other forms of the defor- 
mation energy. Likewise the Morse potential energy, the stacking interaction 
energy experiences also losses which manifests that the stretched H-bonds relax 
and come close to their equilibrium lengths. Subsequent to a drop, the stacking 
interaction energy performs oscillations of fairly large amplitudes throughout 
the time around the mean value of 0.11 eV. Hence, the reduced deformation 
energy of the H-bonds flows into the covalent bond energy and the axial dis- 
tance change energy exhibit analogous temporal behavior, i.e. both rise on the 
average during the redistribution time and approach their quasi-equilibrium val- 
ues. Eventually, the energies E c and Ei oscillate around their respective mean 
value. In other words, the H-bond deformation energy is reduced and its loss 
is deposited into such forms of the energy which contribute mainly to enhanced 
deformations of the covalent bonds. However, the energy associated with the 
axial distance changes remains still very small compared to all other deforma- 
tion energies. We remark that the mean value of the covalent bond energy has 
grown up to 0.06 eV and the peaks of the oscillations are even in the range of the 
Morse potential and the stacking interaction energy giving evidence for major 
distortions of the covalent bonds. 

The energy redistribution processes in the twisted case are illustrated in 
Fig. 6. Compared to the behavior in the preceding ladder case this time the 
Morse potential energy descends already within 2000 ns down to its mean value 
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pointing to a faster equilibration process. Nevertheless, there appear occasional 
bursts of the amplitude. The stacking interaction energy has rapidly dropped 
down from the initial value ~ 0.17 eV to ~ 0.02 eV after only ~ 30 ns. Sub- 
sequently, the stacking interaction energy exhibits oscillations around a mean 
value of 0.55 eV. However, as the covalent bond energy is concerned we find 
that it increases substantially and the eventually attained mean value exceeds 
even that of the Morse potential energy This behavior has to be distinguished 
from the redistribution process in the ladder case where, in the quasi-equilibrium 
regime, the deformation energy is primarily stored as Morse potential and stack- 
ing interaction energy implying that H-bridge deformations play a larger role 
than any other DNA deformations. Moreover, for the twisted case we observe 
also a much more drastic increase in the axial distance change energy and the 
covalent bond energy relative to the ladder case which is associated with inten- 
sified deformations of the helix backbone. Hence, the process of energy redis- 
tribution connected with the approach of a (quasi) equilibrium state proceeds 
differently in the ladder case and the twisted one causing also in the outcome 
distinct structural modifications. 



13 



3.3 Structural changes 

The dynamics of the structural changes in the course of the relaxation process 
is displayed in Fig. 7 and Fig. 8 for the ladder and twisted system, respectively. 
More precisely, we depict the time evolution of the mean value of the lateral 
(radial) deformations 

1 N 

n=l 

the angular deformations 

1 N 

and the mean (axial) distance between two adjacent bases on a strand is given 

by 

1 N 

i=l,2 n=l 

respectively, d (I) determines also the average length change of the hydrogen 
(covalent) bonds. We have directly superposed the graphs from the pulling 
(with an inversed time) and relaxation process. With perspective to the mean 
lateral deformation d we observe that the relaxation process makes fast progress 
relative to the pulling operation so that the mean lateral deformation subsides 
rapidly. Conversely, the mean axial deformation I shows immediately sizeable 
growth in amplitude compared to the corresponding small changes during the 
pulling process. Finally, for the twisted system the mean angular deformations 
do not return to its starting zero value but remain close to the level obtained 
at the end of the pulling procedure such that the helix possesses reduced twist 
relative to the equilibrium configuration. 

Comparing the mean lateral deformations of ladder and the twisted system 
(starting in both cases from the value J(0) = 0.595 A) we conclude that for 
the twisted system d(t) arrives by far much faster at the equilibrium mean 
value d <~ 0.075 A than it is the case for the ladder system attaining the 
mean value d ~ 0.125 A. There remain permanent lateral extensions of the 
ladder and the twisted system deviating from the original equilibrium shape. 
On the other hand, the time evolution of the corresponding variance, Ad = 
l/7V^^ =1 (d„ — J) 2 , reveals that the average deviation from the mean value is 
markedly weaker in the twisted case (Ad < 0.032 A) compared to the consid- 
erably strong deviations in the ladder case (Ac? < 0.084 A). This shows that 
the attractive stacking interaction between successive base pairs gets stronger 
in the twisted system favoring the corresponding inter-base distances to retain 
faster a similar value. In fact, those base pairs which have been forcefully 
separated by the stretching recombine within less than 1000 ns reversing the 
strand separation. Interestingly, the recombination of the DNA duplex (associ- 
ated mainly with the relaxation of the lateral components) takes shorter times 
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Figure 7: Relaxation dynamics: Time-evolution of the structural changes of the 
ladder system, (a) Mean lateral deformations d(t). (b) Mean changes of the 
distances between consecutive base planes l(t). 



than the forced separation. However, the relaxation is not complete because 
the relaxed configuration with its fluctuating positions of the bases still differs 
from the original equilibrium configuration with resting base positions. Further- 
more, distinct from the non-fluctuating behavior during the pulling process the 
relaxation curves exhibit fairly strong fluctuations. Evidently, the relaxation 
(closing) does not take place as the reversed process of the pulling (opening) . 

Analogous to the temporal behavior of the lateral extensions the mean dis- 
tance between two successive bases on a strand l(t) exhibits irregular oscillations 
such that for the twisted system the large-amplitude fluctuations are up to two 
times stronger than the ones for the ladder system. In contrast to the perma- 
nently deformed lateral components (with non-zero mean value d) the average 
axial distance changes fluctuate around the original equilibrium values l n ,i — 0. 
Finally, for the twisted case the mean twist angle performs small-amplitude 
oscillations leaving the helix in slightly distorted angular shape. 

3.4 Comparison between the pulling and relaxation pro- 
cesses 

In summary, the ladder system and the twisted one have in common that when 
both start from an initial partially opened duplex corresponding to a lateral 
elongation pattern they reach due to the relaxation process a quasi-equilibrium 
regime not far from the initial equilibrium configuration. However, unlike for the 
original equilibrium configuration in the newly adopted quasi-equilibrium regime 
the structural coordinates are not at rest but perform oscillatory motions with a 
mean period duration of 50 ns. Correspondingly, this quasi-equilibrium regime 
is characterized by energy exchanges between the various degrees of freedom 
related with the alterations of the H- bridges, the covalent bonds, the longitudi- 
nal elasticity potentials and the stacking interactions. Moreover, the adaptation 
process towards the quasi-equilibrium takes much longer time in the ladder case 
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Figure 8: Relaxation dynamics: Time-evolution of the structural changes of the 
twisted system. Superimposed is the temporal behavior of the deformations 
during the pulling process with an inversed time, (a) Mean lateral deformations 
d(t). (b) Mean changes of the distances between consecutive base planes l(t). 
(c) Mean angular deformation (j>(t). 



than in the twisted one. Hence the ladder system is more rigid with respect 
to spatial deformations than its twisted counterpart which with its additional 
angular degrees of freedom seems to possess more flexibility but less adapt- 
ability. In particular we observe that, with regard to the lateral deformations 
the twisted double-stranded system almost reverts to the original non-distorted 
equilibrium shape on the expense of substantial fluctuations in axial direction. 
In contrast, the ladder system is eventually left in comparatively strong lat- 
eral deformation but without significant axial distance fluctuations pointing to 
strong backbone rigidity. These results demonstrate that the relaxation process 
does not take place as the complete reversal of the pulling process. Similar 
behavior has been observed for the force-induced melting of the DNA double 
helix and its subsequent reannealing 3 . 

Finally we remark that the velocity of the complete recombination of the 
DNA molecule is lower than the velocity imposed by the forcing unit not un- 
common for the opening-closing cycle of DNA molecules (see e.g. [8*]. [24] ) . 
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3.5 Guided relaxation 



We have also studied the elastic response of the stretched DNA double chains 
for the twisted system when at the end of the pulling process the molecule is not 
immediately exempted from the forcing unit. The starting point is an unopened 
molecule whose terminal end is pulled laterally with a velocity v p — 20 fims^ 1 . 
When a certain lateral extension of the terminal base pair is achieved the di- 
rection of the pulling motion is reversed (without waiting time) which returns 
the forcing unit to the position it had before opening allowing for guided relax- 
ation. We monitored the mean lateral extension d(t) for different values of the 
velocity of the backward motion, namely v p = 20/itos -1 , v p = 10/zms -1 and 
v p = 5 pirns' 1 . The results regarding the structural changes are shown in Fig. 9 
for backward velocity Vb = 20 /ims -1 . As expected, at the end of the backward 
guidance of the terminal base pair the whole of the double helix restores virtu- 
ally the original (undistorted) shape. Analogously, the associated deformation 
energies decrease steadily in the course of the relaxation process and adopt their 
initial zero content which they had before the pulling operation. However, we 
observe an asymmetry between the forward and backward elongation curves. 
Particularly, in an early stage of the retraction of the pulling unit the corre- 
sponding relaxation curves are in advance of their counterparts for the pulling 
process which seems to present a typical behavior of DNA molecules (see e.g. 

0,0, na) 

With lowered backward velocity the process till complete relaxation is achieved 
takes correspondingly longer times. 

3.6 Parameter dependence 

For a parameter study we have varied the stacking interaction strength 5* and 
simulated the pulling process to follow afterwards the relaxation dynamics. As 
the main result we found that the larger S is the more rapid proceeds the 
relaxation process, viz. the faster is the quasi-equilibrium regime attained. 
Furthermore, it holds in general that the energy migration takes more time in 
the case of the ladder system than for its twisted counterpart. 

With view to a dependence on the molecule length N, i.e. the number of 
base pairs, we remark that the time scale for the relaxation process remains 
basically the same when N is varied in the range from 10 to 30 providing strong 
evidence for the enormous elasticity of DNA molecules on the lower amplitude 
scale. In contrast, for large-amplitude elongations the dissociation time of DNA 
molecules increases from a fraction of a second to hours over the range of N = 10 
to 20 base pairs [T2"] . 

In addition, with regard to asymmetry we investigated also cases when the 
ends of the strands are pulled laterally with unequal forces (in the extreme 
case only one strand is pulled) resulting in elongation patterns being asymmet- 
ric with respect to the central helix axis. We observed then that the systems 
behave during the subsequent relaxation dynamics like in the symmetric cases 
discussed above and attain finally quasi-equilibrium regimes being not far away 



17 




1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 



Figure 9: The relaxation dynamics for the twisted system for backward velocity 
Vb = 20 [irajs. (a) Mean lateral deformations d(t). (b) Mean changes of the 
distances between consecutive base planes l(t). (c) Mean angular deformation 

m 



from the respective original equilibrium configurations with fluctuating struc- 
tural coordinates. 

4 Summary 

In this paper we have demonstrated that within the frame of an oscillator net- 
work model we are able to study numerically the mechanical stability and elas- 
ticity properties of DNA molecules. Our model of coupled oscillators takes 
essential microscopic degrees of freedom of DNA and the inherent interactions 
between them into account. We have focused our interest on the opening-closing 
dynamics of double-stranded DNA molecules. The opening of a DNA molecule 
has been forced by mechanical stress imposed at a terminal end of the molecule 
bringing it into a non-equilibrium state. We have followed the subsequent relax- 
ation towards an equilibrium state with recombination to its double-stranded 
conformation. The similarities and differences between the relaxation dynamics 
for a planar ladder-like DNA molecule and a twisted one have been discussed. 
In particular we have shown that the attainment of a quasi-equilibrium regime 
proceeds faster in the case of the twisted DNA form than for its ladder coun- 
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terpart. This has led us to the conclusion that in the twisted form the DNA 
molecule is more flexible with respect to conformational changes than in the 
planar ladder-like version. 

There are some limitations in our model approach and for an improvement 
one should certainly consider the dependence on the DNA sequence, the impact 
of the chemical environment and the influence of temperature on the mechanical 
stability. 
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